Origin of giant bulk Rashba splitting: Application to BiTel 
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We theoretically propose the necessary conditions for realization of giant Rashba splitting in 
bulk systems. In addition to (i) the large atomic spin-orbit interaction in an inversion-asymmetric 
system, the following two conditions are further required; (ii) a narrow band gap, and (iii) the 
presence of top valence and bottom conduction bands of symmetrically the same character. As a 
representative example, using the first principles calculations, the recently discovered giant bulk 
Rashba splitting system BiTel is shown to fully fulfill all these three conditions. Of particular 
importance, by predicting the correct crystal structure of BiTel, different from what has been 
believed thus far, the third criterion is demonstrated to be met by a negative crystal field splitting 
of the top valence bands. 
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Exploiting the spin degree of freedom of the electrons 
is one of the primary goals in the rapidly growing field 
of spintronics. A promising candidate to achieve this 
goal is so-called Rashba effect, which relies on the spin- 
orbit interaction (SOI) of the carriers in an inversion (/) 
asymmetric environment |^ . This effect has been exper- 
imentally observed for a number of non-magnetic metal- 
lic surfaces [^-Q and also demonstrated to exist at the 
interface of the semiconductor hetrostructures For 
most of these systems, the level of Rashba spin split- 
ting (RSS) are found to be rather small (at most several 
meV). However, there have been a few exceptions, e.g. 
the Bi-covered Ag(lll) surface, for which the angle re- 
solved photoemisson spectroscopy (ARPES) has revealed 
a giant RSS of the order of 200 meV Following these 
discoveries, attempts have been made to realize RSS in 
three-dimensional systems, as such systems are expected 
to be an ideal laboratory for exploring many novel phe- 
nomena, e.g. the spin Hall effect and resonance enhanced 
magneto-optical conductivity. 

BiTel, a polar layered semiconductor, has been very 
recently revisited from this point of view. The ARPES 
measurements clearly show a gigantic RSS among the 
lowest conduction bands (LCB's) in the bulk BiTel 
leading to a substantial shift, kcBM = ±0.05 A-i, in the 
position of conduction band minimums (CBM's). The 
level of spin splitting obtained at CBM astonishingly 
reach to ^ 0.4 eV lying among the highest discovered so 
far. The corresponding Rashba energy Eji — Eqbm^Eq, 
where Eq indicates the energy of the two LCB's at their 
crossing point, is found to be over 100 meV. The spin- 
resolved ARPES measurements further reveals that these 
bands are fully spin-polarized, reassuring the entire sys- 
tem is indeed subject to a giant RSS 

Motivated by this discovery, we have performed a the- 
oretical study based on the perturbative k • p formalism 
and backed by the first-principles calculations and group 
theoretic analysis to investigate the origin of giant RSS 



in bulk materials. Through this study, we have identified 
three conditions required for realization of this intriguing 
phenomenon. These conditions are shown to be closely 
related to the relative ordering and symmetry character 
of the bands near Fermi level, Ep. As a representative 
case, BiTel is shown to fully meet all these conditions, 
owing to its unusual electronic band structure near Ep. 

As a starting point, we describe the general k-p Hamil- 
tonian via perturbation theory (PT). Given the solution 
-ff(ko) at k = ko, it can be expressed for nearby k as. 
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Here, V , a and p denote the crystal potential, Pauli 
matrices and the momentum operator, respectively, and 
q = k — kg. Considering only the linear-in-A: spin split- 
tings, one can show that they can arise due to iJ^^^ by 
the use of the first-order PT or the coupling between per- 
turbing terms ;^q • P and H^'^^ in the second order PT. 
The spin splitting arising from i/^^^ is, however, expected 
to be much less than that coming from i?^^-', and hence, 
unlikely to cause any giant RSS (i.e. > 100 meV). This 
is because H^^^ {H^"^^), as refiected by its k-dependence 
(p-dependence), originates from the crystal (atomic or- 
bital) momentum. Since the velocity of the electron in 
its atomic orbit is far greater than the velocity of a wave 
packet, the spin splitting is accordingly expected to be 
strongly dominated by if^^^ 0. The respective second 
order perturbative correction in energy is given by ^ , 
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FIG. 1. (Color online) (a) Crystal structure of BiTel and (b) 
simulated XRD pattern for the optimized and non-optimized 
structures. 



where Ui and Si are the eigenstate and eigenenergy corre- 
sponding to the state i at kg, respectively, and cc. stand 
for the complex conjugation. 

Equation (^) clearly indicates that the level of spin 
splitting is directly dependent on three conditions: (i) 
the strength of the spin orbit interaction (represented by 
if^^-') (ii) the energy difference between the neighboring 
states m and n, and (iii) the symmetry character of their 
corresponding eigenstates, determining if (itml^f'-^-' |uti) 
is symmetrically allowed or not. In brief, states ener- 
getically close to each other and symmetrically of the 
same character, can effectively couple with each other, 
and hence produce a large spin splitting, of course if the 
host atoms maintain a strong SOI. The first condition 
is of course a rather obvious requirement, already well- 
known in the context of surface RSS . The other two 
are, however, less trivial and require more attention. In 
the case of semiconducting bulk materials, they can be 
satisfied if the band gap is sufficiently narrow and, more 
importantly, if both LCB's and the highest valence bands 
(HVB's) are symmetrically the same. The last condition 
requires an anomalous ordering of the bands near E p 
which is not usually allowed in conventional semiconduc- 
tors. However, some polar semiconductors can excep- 
tionally meet this criterion due to the negative crystal 
field splitting (CPS) of their top valence bands (TVB's). 
BiTel is one such example, which in the following will 
be shown to fulfill all these three conditions and thereby 
exhibiting a giant bulk RSS. 

Having a trigonal structure with the space group P3ml 
(No. 156), BiTel is a polar compound in which Bi, Te 
and I form stacking layers along the c-axis. A charac- 
teristic feature of P3ml group, as shown in Fig. ||(a), is 
the lack of /-symmetry. The highest symmetry opera- 
tion allowed is C^v along the c axis. The corresponding 
experimental lattice parameters a and c are 4.339 A and 
6.854 A, respectively |jl^. Assuming Bi to be at ori- 
gin, Te and I have been experimentally proposed to be 
respectively at (2/3,1/3,0.6928) and (1/3,2/3,0.2510) 
sites [|ll|. However, after a full structural optimization 
of atomic positions [n2|, we have surprisingly found that 



Te and I change their positions to (2/3, 1/3,0.7482) and 
(1/3,2/3,0.3076), respectively. In other words, our pre- 
dicted Bi-Te (Bi-I) distance is exactly equal with the ex- 
perimentally proposed Bi-I (Bi-Te) distance. Such a cor- 
rection, as will be discussed later, leads to a number of 
fundamental changes in the electronic structure of BiTel, 
including the appearance of a gigantic RSS, similar to 
what was observed by ARPES [§. The failure of XRD 
experiment might be due to the fact that Te and I both 
have nearly the same ionic radii (133 pm and 131 pm, re- 
spectively ) and atomic charges (52 and 53, respectively). 
Accordingly they likely produce rather indistinguishable 
features in the XRD pattern. Indeed the simulated XRD 
patterns [|l3| , as shown in Fig. 0(b) , turn out to be nearly 
identical for the both structures. Thus, the utilization of 
more sophisticated experimental techniques seems to be 
necessary for the proper identification of atomic positions 
in BiTel. In the rest of this letter, we compare in detail 
the electronic structures of both the non-optimized and 
optimized BiTel , to asses if the above mentioned condi- 
tions are indeed required for a giant bulk RSS. 

Figure || shows the respective band structures for the 
both systems. In the absence of SOI, they show a semi- 
conducting behavior with an energy gap Eq ^ 1.2 eV. 
As shown in Fig. ||(a) and (b), the lowest Eq is com- 
monly found to be not at the BZ center, but at A point 
where — ky — and fc^ = tt (see Fig. ^(f)). Up 
to 4 eV above Ep, the (six) conduction bands are pre- 
dominated by Bi-6p states, whereas the Te-5p and l-5p 
most strongly contribute to the (twelve) valence bands 
down to -5 eV below Ep- Another important similar- 
ity is the presence of a rather large CFS among these 
bands. Without SOI, for any fc-point other than those 
along T-A or the ones with accidental symmetry, all the 
p-type bands split into doubly degenerate bands. Along 
T-A, such states are allowed to form either 2-fold or 4- 
fold degenerate bands. Such a trend of CFS can be well 
described using the group theory. Within P3to1 group, 
at the BZ center and along F-A all the fc-points have C^y 
symmetry. Without spin, all the bands at A are thus 
transformed according to one of the single group repre- 
sentations of Csi,, (a complete character table for C^^ 
can be found in the supplemental file). Of our particular 
interest are {s,pz} — Ai and {px,Py} — ^3, explaining 
why the p-type valence and conduction bands along T-A 
are either two-fold or four-fold degenerate. For the other 
internal BZ points all the single representations are one 
dimensional and hence non-degenerate except for spin. 

Comparing the ordering of CFS of conduction bands 
at A, one can notice a similar trend, that is, A1-A3 in 
the increasing order of energy. The same trend holds for 
TVB's in the non-optimized structure. However, in the 
optimized structure the ordering of TVB's is opposite, i.e. 
^3-^1 (see Fig. ^(b)). This accordingly implies the exis- 
tence of a negative CFS among this group of bands. As a 
result, both the LCB's and HVB's turn out to be Ai (p^) 
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FIG. 2. (Color online) Calculated electronic band structures for the non-optimized structure (a) without SOI and (c) with 
SOI and for the optimized structure (b) without SOI and (d) with SOI. (e) An scaled-up view of band dispersions along H-A-L 
direction for the optimized structure, (f) Corresponding high symmetry fc-points in the hexagonal Brillouin zone of BiTel. In 
(a) and (b), the oblique numbers indicate the band degeneracies at A point for each structure (see the related discussion). 



type and 2-fold degenerate. It is to be noted that such a 
negative CFS has already been observed for a number of 
chalcopyrite-type semiconductors, e.g. CdSnP2 [Q and 
CuAlS2 ^ as well as AIN [|6). The reason has been 
attributed to the strong ionicity of the atomic bondings, 
leading to a substantial structural distortion along their 
high symmetry axis. 

Turning on SOI, the band structures undergo yet an- 
other drastic change. As shown in Fig. ^(c) and (d). 
Eg is closed down to ^ 0.28 eV, nearly 5 times smaller 
than that obtained without SOI. This is mainly due to 
the strong SOI of Bi which shifts downward the j — 1/2 
bands by nearly —2Aso, where Ago denotes the atomic 
SOI energy (for Bi it is found to be ~ 0.5 eV). De- 
spite this similarity, the trend of spin splitting near Ep 
is strikingly different between the two system. A com- 
parison between Fig. ||(c) and (d) clearly reveals that in 
the optimized structure a huge RSS takes place among 
both LCB's and HVB's at A point, whereas the non- 
optimized structure fails to yield such a feature. As 
shown in Fig. ^(e), in the optimized structure, CBM and 
VBM are both shifted by nearly kcBM — kvBM ~ ±0.05 
from A in the {kx, ky) plane with an Ej^ = 113 meV, 
in excellent agreement with the ARPES data ^] . Such a 
good agreement is a strong indication that our predicted 
structure is indeed correct. 

Having confirmed the existence of a giant RSS in bulk 
BiTel, we next address the main question: whether this 
effect arises due to the fulfillment of all the three condi- 
tions, pointed out earlier. As already might be under- 
stood, the both optimized and non-optimized structures 
meet the first two criteria, namely the strong SOI in an 
/-asymmetric environment and narrow band gap. How- 
ever, the last condition is only satisfied in the optimized 
structure as its both LCB's and HVB's are symmetri- 
cally the same. To be more specific, upon introduction 



of SOI, the previously defined single group representa- 
tions Ai and A^ transform to their double group coun- 
terparts (see the character table in supplemental file) . In 
Cat, space the transformation is such that Ai — )■ and 
A3 — )■ ^40^50^6. In other words, Ai transforms to A4, 
whereas A3 is split into two 2-fold bands A4 and A.5 © Aq 
( the latter is hereafter simplified as A^^). Since TVB's 
of optimized structure undergo a negative CFS, the cor- 
responding HVB's and LCB's are both of A^ character. 
Fig. H schematically shows the effects of CFS and SOI. 

Turning back to Eqs. (|l|) and (^, one can use the 
method of invariants to find an effective spin-splitting 
Hamiltonian linear in terms of k. For C3t, symme- 
try, it is easy to show that, the only possible choice is 
Hq oc {(Txqy — <Jy(lx) bccausc both {qx,qy) and {ax,cry) 
belong to the A3 representation. Hq is clearly a Rashba- 
type Hamiltonian, giving the following spin-split energies 

Aem — ±am\/qx + qy- ^s described in the beginning. 
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FIG. 3. (Color online) Diagrammatic representation of band 
splitting due to the crystal field splitting (CFS) and spin-orbit 
interaction (SOI) and their combination in (a) non-optimized 
and (b) optimized BiTel. 
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thus pursue with PT to quahtatively determine the sec- 
ond order correction in Rashba parameter am, denoted 
as am ■ With the help of group theory, it turns out 
that in Eq. (|) {um\H(-'^'>\un) 7^ 0, if and only if both 
Um and Un are of ^4 character, implying that linear- 
in-k splitting is symmetrically forbidden for A^^q bands. 
Accordingly, we just need to consider Ui = and A41 
(see Fig. ^). To derive an analytical form for am"^ from 
Eq. (H), analogous to Gutche and Jane basis ^ 0, 

we define IA4 U) = ±\/^^\{x±iy) ;t)+94|^tl), and 

\Ai, U) - ±^\{x±iy) ;t> - V^~gI\Z U) , where 

is to enforce the orthogonality of the basis sets and + (— ) 

sign corresponds to spin f (spin J,). The Ai-like basis \Z) 

is defined as \Z) — qs\s) + qz\z) with _|_ ^2 _ ^ 

To simplify our derivations, we recall the fact that the 

(2) 

states with small {Sm — En) dominate AeJn (k) and, thus 
am\ From Fig. ^(c) and (d), it is evident that |e4c — 
£4'c| > 1.30 eV and \e4v — £4't,| > 1.15 eV. On the other 
hand, for the optimized and non-optimized structures, 
the respective Ae4_4 — e^c—^Av and Ae4^4/ = e4c— e4'„ are 
both below 0.5 eV and very close to their Eq. Thus, as a 
good approximation for the former (latter), a^^^ = a^^J^^ 

{af^ = aic!4't,) with, 



aicUv = ^^^^^ [g4t;y 1 - qlcPx,Z + '74c Y 1 - <llvPz,x] 

where, A,j = {A,c t W^^^lA^y f) and P^^.z = 
—ih{xc\d/dx\Zy). 

Constructing a set of maximally localized Wannier 
functions [p^-pO| for the conduction and valence bands 
around Ep, we have estimated g|„ and ql^ for the op- 
timized and non-optimized structures. For the latter, 
ql^ — 0.55 and q^^ — 0.43, both very close to the critical 
value 0.5. Assuming Px,z — Pz,x, one can immediately 
find from Eq. (||) that a^^^ ~ 0. Since a^?^ 4^ = -o^^aIw^ 
then a^'v is also expected to be nearly zero. This clearly 
explains why the non-optimized structure shows almost 
no spin splitting among its HVB's and LCB's near the 
A point. As for the optimized structure, the situation is 
completely different. Here, q\^ — 0.886 and q\^ = 0.5, im- 
plying that the HVB's are predominantly .Z-type. Con- 
(2) (2) 

sequently, a\^ and a\^ turn out to have appreciable val- 
ues with equal magnitudes but opposite signs. In other 
words under this situation, for both HVB's and LCB's, 
the absolute value of a can be nearly the same, but their 
signs are always opposite. That's exactly what we can 
see in the optimized BiTel as it also shows similar trend 
of spin-splitting among its LCB's and HVB's such that 
their corresponding kcBM and kvBM are almost at the 
same place. Here, It is important to emphasize that such 
a second order perturbative RSS is a direct result of an 



(i) anomalous ordering of top valence bands due to the 
existence of negative CFS which allows the adjacent ^4^ 
and A4C to be symmetrically of the same character and, 
hence, to be coupled with each other thorough a per- 
turbative Rashba-like Hamiltonian. Due to (ii) the large 
SOI of Bi leading to (iii) substantial band gap narrowing 
, such a coupling can be effectively very strong. These are 
the three key factors for realization of giant RSS in BiTel, 
and very likely any other giant bulk Rashba splitting ma- 
terial. For the other candidates a negative CFS can either 
intrinsically exist due to the strong anisotropic ionicity 
of their atomic bondings or be externally produced e.g. 
through a pressure-induced structural distortion. 

In summary, in this study we combined the first- 
principles calculations with a group theoretic analysis to 
investigate the origin of giant Rashba splitting in bulk 
systems. As a representative case, It was shown that in 
BiTel the interplay between the giant SOI of Bi and effec- 
tively large negative CFS of the TVB's led to a substan- 
tially strong coupling between the narrowly separated 
HVB's and LCB's via a perturbative Rashba-like Hamil- 
tonian. Such conditions were expected to be vital for 
realization of giant RSS in other bulk candidates, also. 
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ogy (FIRST program) on "Quantum Science on Strong 
Correlation" . 



* bahramy@riken.jp 
[1] E. I. Rashba, Sov. Phys. Slid State 2 , 1109 (1960). 
[2] S. LaShell et al, Phys. Rev. Lett. 77, 3419 (1996). 
[3] C. R. Ast et al., Phys. Rev. Lett. 98, 186807 (2007). 
[4] Y. M. Koroteev et al, Phys. Rev. Lett. 93, 046403 
(2004). 

[5] J. Nitta et al, Phys. Rev. Lett. 78, 1335 (1997). 
[6] K. Ishizaka et al., Unpublished. 
[7] E. O. Kane, J. Phys. Chem. Sohds 1, 82 (1956). 
[8] L. C. Voon, et al, Phys. Rev. B 53, 10703 (1996). 
[9] L. C. Voon, and M. Willatzen, The k.p Method, (Springer, 
Berlin 2009). 

[10] M. Nagano et al, J. Phys. Condens. Matter 21, 064239 
(2009). 

[11] A. V. Shevelkov et al, J. Solid State Chem. 114, 397 
(1995). 

[12] All the calculations, reported here, are carried out using 
PBE exchange-correlation functional and the augmented 
plane wave plus local orbitals method with inclusion of 
SOI as implemented in wien2k program (P. Blaha et 
al, URL: http://www.wien2k.at). The muffin tin radii 
are set to Rmt = 2.5 bohr for all the atoms and the 
maximum modulus for the reciprocal vectors Kmax is 
chosen such that RhixKmax ~ 7.0. Accurate BZ sam- 
plings are performed using a 20 x 20 x 20 k-mesh. For the 
structural optimization, the stopping criterion is when 
the magnitude of force on all ionic sties is less than 0.1 
mRy/bohr. 

[13] Simulated by CRYSTALDIFRACT program, Ver. 5.2.0: 



5 



http:/ /www.crystalmaker.com/crystaldifFract/index.html 
[14] J. L. Shay et al, Phys. Rev. Lett. 24, 1301 (1970). 
[15] V. .Jayalakshmi et al, J. Phys. Chem. Sohds 67, 669 

(2006). 

[16] Y. Tanoyasu et al, Appl. Phys. Lett. 90 261911 (2007). 
[17] E. Gutsche, and E. Jahne, Phys. Status SoUdi 19, 823 



(1967). 

[18] L Souza et al, Phys. Rev. B 65, 035109 (2001). 
[19] A. A. mostofi et al, Comput. Phys. Commun. 78, 685 
(2008). 

[20] J. Kunes et al, Comput. Phys. Commun. 181, 1888 
(2010). 



